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Metallic split-ring resonators possess dominant electric dipolcs as well as considerable magnetic dipoles under 
proper excitations. Full-wave numerical approaches are frequently employed to simulate adjacent split-ring 
resonators, but simulations cannot explain the underlying physics. An analytical theory based on a dynamic 
electric dipole approximation is developed here. Detailed theory-simulation comparisons demonstrate that 
this theory can quantitatively describe the interaction strength of coupled split-ring-resonators under certain 
circumstances. © 2010 Optical Society of America 
OCIS codes: 220.1080, 240.6680, 350.4990 
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Since its introduction in 1999, metallic split-ring res- 
onators (SRRs) have attracted intense attention because 
of their important roles in artificial metamaterials [1]. 
Under appropriate illumination, circulating current will 
be excited inside the SRR and lead to magnetic dipoles 
comparable in magnitude to electric dipoles, making 
metallic SRRs ideal building blocks to achieve strong 
magnetic responses in the optical region. A single SRR 
consequently is referred to as a magnetic photonic atom 
[2]. Recently, several experiments are devoted to mul- 
tipole coupling between SRRs, especially electric and 
magnetic dipole interaction [3-6]. For instance, planar 
metallic SRR arrays are fabricated with different lat- 
tice spacing along different directions, and both mag- 
netic and electric near-field dipole coupling are found 
to simultaneously influence the spectral positions of the 
plasmonic resonances [5]. Another experiment studied 
isolated metallic SRR dimers with different configura- 
tions. Changing the relative magnetic dipole phase of 
the constituent particles identifies the strength of mag- 
netic dipole coupling [4]. Two methods are frequently 
employed to interpret the experimental observations, 
full- wave numerical simulations and quasistatic coupling 
models. Simulations are able to reproduce the experi- 
mental measurements quantitatively, but shed little light 
on the underlying physics [4]. The quasistatic model, on 
the other hand, obtains variables by fitting experimental 
data. Consequently it can not make an a priori predic- 
tion [5-7]. 

In this Letter we will show that a priori predictions of 
the electromagnetic interaction are feasible for specific 
SRRs structures. More particularly, we carefully com- 
pare the dynamic radiation emitted from the magnetic 
dipole of a single SRR with that from its electric counter- 
part. It is found that at particular "safe" spatial regions 
the electric dipole exclusively dominates the radiation 
field. An individual SRR consequently can be approxi- 
mated by an isotropic electric dipole whose polarizability 
is determined by the particle geometry as well as the il- 
lumination. This approximation can be further extended 



to multiple SRRs where each constituent particle sits in 
the safe zone of other particles. Two dimers are studied 
to verify our theory. By comparing the analytical results 
with the corresponding full- wave numerical results, it is 
demonstrated that our theory allow us to predict the 
coupling strength quantitatively. 

We start with the optical response of a gold SRR as 
shown in Fig. 1(a), which closely matches the experimen- 
tal samples of Ref. [5]. It is excited by an x-polarized 
plane wave propagating in the z direction. The absolute 
extinction and absorption cross section of the particle 
are shown in Fig. 1(b) [8]. A peak corresponding to the 
fundamental plasmonic resonance was found at a wave- 
length of 1352 nm. Because of the structural symme- 
try breaking, a z-component magnetic dipole m^ is gen- 
erated together with a dominant x-component electric 
dipole dx- Numerical simulations further suggest that 
the ratio /3 = —mz/dx/c is around 0.29? [8]. The mag- 
netic dipole therefore has an amplitude comparable to 
its electric counterpart. In addition, because the extinc- 
tion spectra are very sensitive to the illumination polar- 
ization, both electric and magnetic dipolcs arc strongly 
anisotropic [5,8]. 

It is well known that the electric field, in the near and 
far zone, for an electric dipole source p is given by [9] 



Ec(r) = [A{kr)l + B(fcr)T] • p 



(1) 



with k being the wave number, I being an unit tensor of 
rank 2. The tensor T has component as Tij = riTj/r"^, 
and the coefficients are written as [10] 
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The electric fields emitted from the electric dipole d^ 
of a single SRR are therefore direction-dependent: The 
A term survives for any radiation direction while the B 
term vanishes when r parallels the y axis. In a similar 
way, the radiation for a magnetic dipole source m bears 
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Fig. 1. (a) The x-coniponent of the current distribution 
of an individual spht-ring resonator at its fundamental 
resonance. The incident pulse is x-polarized and propa- 
gating in the z direction, (b) The extinction and absorp- 
tion cross sections of the single split-ring resonator. The 
frequency is normalized to the fundamental resonant fre- 
quency Wo- The bulk plasma frequency of gold is taken 



as Wp = 1.367 X 10^'' 
frequency 7 = 6.478 x 10 



the phenomenological collision 
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the form 
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with Zq = \/iM)/eo being the vacuum impedance. The 
polarization of the radiation field therefore is perpendic- 
ular to the plane defined by r and m. Consequently, the 
electric field emitted from the magnetic dipole niz is po- 
larized in the xy plane and is perpendicular to r when r 
lies in the xy plane. 

We now can analyze the radiation field of the single 
SRR in the xy plane. Because its fundamental wave- 
length, 1352 nm, is much larger than its characteris- 
tic size, higher-order multipoles fall off rapidly and the 
radiation emitted from the structure will come mainly 
from its electric and magnetic dipoles when kr > 1 [9]. 
More specifically, when r parallels e^, Eo is x-polarized 
and E,n is y-polarized, and the ratio of their amplitude 
Ec/Eni equals 2/|/3|fcr. The magnetic dipole then domi- 
nates the radiation zone where fcr ^ 1, and the electric 
dipole dominates a limited region where 1 < kr < 1/|/3|. 
Similarly, the polarization of both Ee and E^ are along 
the X axis for r parallel to Sy. The amplitude ratio 
Ec/Et, with Et being Ec -I- Em, is further given by 
[kr + i{l - fc^r^)] / [kr + \13\ + i{l - fc^r^)] . As a resuh, 
we can safely neglect the magnetic dipole and use the 
dynamic electric dipole approximation when kr ^ |/3|. 
Similar conclusions apply for r along the z axis. The 
electric vector Em is exactly zero because r parallels m^ , 
and the Eo is x-polarized. For a considerably large fcr, 
the magnetic dipole has no infiuence on the radiation 
field. 

The discussions above can be extended to multiple 
SRRs, as long as each constituent particle sits in the 
public "safe" zone. To validate the theory, we investigate 
two different configurations of SRR dimers, side-by-side 



and on-top, which have been studied experimentally [4]. 
The dipole polarizability a(cj) of an individual SRR is 
determined first, which connects to the extinction cross- 
section spectrum cr(cj) as [8] 
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where Eq stands for the incident electric field, J repre- 
sents the polarization current and the integration is per- 
formed over the particle volume. We neglect the retarded 
effect because the particle height, 25 nm, is much shorter 
than the incident wavelength. Moreover, a Lorentz model 
describes the polarizability a{u!) around the resonance 
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where ujq is the resonant frequency. / stands for the os- 
cillator strength and r measures the phenomenological 
damping force [9]. For the SRR shown in Fig. 1(a), we 
find that t/loq « 0.072 as well as f/r w 1.29 x 10"^. 
It should be emphasized that the extinction spectrum is 
experimentally measurable [4,11] and the parameters of 
the Lorentz model can then be obtained by fitting the 
experimental results. 

By approximating each particle of an SRR dimer as 
an electric dipole, a set of coupled equations describes 
its optical properties 

Pi = a{iu) [Eo + Ec(d, P2)] , P2 = a{Lo) [Eq + Ec(d, pi)] , 

(6) 
where d stands for the distance between pi and P2, with 
Pi and P2 being the dipole moment of each particle, re- 
spectively. The symmetry of the equations immediately 
implies that the incident plane wave Eq only excites 
the symmetrical (in-phase) mode in either configuration. 
The extinction spectrum can be obtained by solving the 
set of equations. For the side-by-side structure it is given 

by 
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Because A and A+B tend to zero when d goes to infinity, 
the equations above lead to an intuitive fact: The total 
power taken from the incident wave by two identical par- 
ticles without coupling is twice of that by a single parti- 
cle. Furthermore, the resonant frequency a;^ corresponds 
to the vanishing of the real part of the denominator of 
CTe [10], we therefore have 



ineo Rc[A{ujrd/c)] = - 
for the on-top configuration, and 

47reo Rc[A{uird/c) + B{uJrd/c)] 



(9) 



(10) 
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Fig. 2. Two different configurations, sidc-by-side and on- 
top, of two split-ring resonators. The resonant frequen- 
cies are normalized to wq, the fundamental resonant fre- 
quency of the individual split-ring resonator. The extinc- 
tion cross sections are normalized to the as{uJo), the 
extinction cross section of the single split-ring resonator 
at its fundamental resonance. The analytical and nu- 
merical results arc shown with solid and dotted curves, 
respectively. 



for the side-by-side configuration. Notice that the damp- 
ing force r of the Lorentz model does not appear since 
it mainly influences the bandwidth of the resonance [9] . 
These equations can be solved easily in the near-field 
zone where kr <^ 1. We find that 47reoRc[^] w — l/d'^ as 
well as 47reoRe[y4 + B] ^i 2/d'^, the resultant ujr is there- 
fore bigger than wq for the on-top dimer while smaller 
than luq for the side-by-side arrangement. In addition, 
the resonance shift jw^— wq] of the side- by-side structure 
is considerably bigger than that of the on-top configura- 
tion [12]. The former structure therefore dominates the 
nearest neighbor coupling of equally oriented SRRs in 
dense square arrays [5]. We want to stress that all the 
predictions here are qualitatively consistent with the ex- 
perimental observations [4,5]. 

Equations (9) and (10) arc solved for the separation d 
ranging from 200 to 500 nm, and the results arc plotted 
in Fig. 2(a) with solid curves. The resonant frequencies 
LOr are normalized to wq of the individual particle. A 
finite-difference time-domain method is then applied to 
numerically simulate the two dimers, and the calculated 
results are plotted in Fig. 2(a) with dotted curves [13]. 
The center-to-center distance d, as marked in Fig. 2, is 
gradually increased from 225 nm to 475 nm (correspond- 
ing to kd from 1.0 to 2.2) with an increment of 25 nm. 
Notice that the width of the gap between these two par- 
ticles is given by d — 200 nm; they arc very close to 
each other. The numerical result of the on-top dimer is 
found to agree perfectly with its analytical counterpart 
when d > 250 nm, while a slight discrepancy appears 
for the side-by-side configuration: The numerical result 
oscillates around its analytical counterpart with an av- 
erage relative difference smaller than 1%. It is possibly 
induced by the non-uniform distribution of the polariza- 
tion current inside the SRR. As shown in Fig. 1(a), the 



current is found to be strongly localized around the two 
inner corners. Assuming the current "hot" points are the 
positions of the electric dipoles, the dipolc-to-dipole dis- 
tance is shortened to roughly d— 100 nm for the side-by- 
side dimer. Higher-order multipoles such as an electric 
quadrupole may interfere with the electric dipole to a 
certain degree. 

Although our theory describes the coupling strength 
quantitatively, it can not be extended to the extinction, 
absorption or scattering spectra because these quanti- 
ties need the information regarding radiation along all 
the directions. However, its predictions can serve as a 
first approximation since the dominant electric dipole is 
contained. Take the two SRR dimers studied above as 
examples: their normalized extinction cross sections ver- 
sus d are computed analytically and numerically, and 
the results arc plotted in Fig. 2(b). Unlike the resonant 
frequency, the analytical results of the extinction are 
considerably different from their numerical counterparts, 
with an average relative difference around 10%. In addi- 
tion, the theory must be employed cautiously to study 
the gap-to-gap and back-to-back configurations consid- 
ered in Ref. [4]. The constituent particle in either dimer 
has in-phase electric dipoles but out-of-phase magnetic 
dipoles; Equation (6) then loses its validity. It is still a 
good approximation when d is big enough because the 
magnetic dipole coupling falls off rapidly with increasing 
intcr-particlc spacing. Finally, note that magnctoinduc- 
tive coupling will appear in the side-by-side dimer when 
d ^ c/\/3\ujQ, as suggested by our theory [2]. 

To summarize, by carefully analyzing the fields emit- 
ted from the electric dipole and magnetic dipole of an 
individual split-ring resonator, we demonstrate that in 
certain circumstances a simple dynamic electric dipole 
can be employed to approximate the split-ring res- 
onator. Two configurations of coupled split-ring res- 
onators are studied to validate our theory by detailed 
theory-simulation comparisons. It is shown that this the- 
ory can predict coupling strength quantitatively and 
other optical quantities such as extinction cross section 
qualitatively. 
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